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We report a study of the B meson decays, _B+ J/i>(j)K+, B" ~+ J/i>(j)Ks, B° J/ip(t>, 
B° — > and B" —^ J/iprj' using 56 million BB events collected at the T(45') resonance with the 

BABAR detector at the PEP-II e'^e~ asymmetric-energy storage ring. We measure the branching 
fractions B{B+ J/'ilj(f)K+) = (4.4 ± lA{stat) ± Q.b{syst)) x lO"'^ and S(B° ^ J/il^(t>Ks) = (5.1 ± 
1.9{stat) ±0.5{syst)) x 10~^, and set upper limits at 90% confidence level for the branching fractions 
B{B° ^ J/iP<l)) < 9.2 X 10-^ B{B° J/iPr;) < 2.7 x 10-^ and B{B° -> J/ipri') < 6.3 x 10-^ 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 

Recent observations of the B meson decays B J/ipir Q and J/ipp are evidence for the Cabibbo- 
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suppressed transition b ccd via the color-suppressed 
diagram shown in Fig. ^ (a). Here we present a search 
for color-suppressed modes with hidden strangeness, ss, 
in the final state: B J/ipV^ J/''Pv' ^ J/'4'4' ^^nd J/tpcf)K. 
The decays -B° — > J/tprj and J/ipv' occur via the 

same diagram, Fig. (a), and should have a rate com- 
parable to B — > J/V'TT. If large enough samples can be 
isolated, these CP eigenstates could be used to test CP 
violation Models based on the heavy quark factor- 
ization approximation by A. Deandrea et al. Q are used 
to predict that the branching fraction for — > J/ip7] 
is a factor of 4 smaller than that for i?" J/t/itt". As- 
suming that the decay B^ J/'^4> is a color-suppressed 
mode with rescattering as shown in Fig. (b) , then the 
absence of a signal would indicate that the rescatter- 
ing effects are negligible. The decay B — > J /ip(f)K is 
a Cabibbo-allowed and color-suppressed decay via the 
transition hq ccsssq, where the ss quark pairs are 
produced from sea quarks or are connected via gluons as 
shown in Figs.n(c) and (d), respectively. This particular 
three-body decay would be of interest in the search for 
hybrid charmonium states that decay to the final state 
J /^4> |5| . In this paper we report on branching fractions 
or upper limits for J/tprj, J/^rj' ^ J/ipff'^^y and 

J/ib<j>Kl 



(a) 



(b) 




(0) 



c 




- J/y 



FIG. 1: Quark diagrams: (a) tree diagram for B J/tpn and 
J/ipp, (b) rescattering for B J/^4>, (c) strange sea quarks 
and (d) gluon coupling for B — > J/tpifiK. 

The data used in this analysis were collected at the 
PEP-II asymmetric-energy e"'"e~ storage ring with the 
BABAR detector, fully described elsewhere [^. The 
BABAR detector contains a five-layer silicon vertex tracker 
(SVT) and a forty-layer drift chamber (DCH) in a 1.5-T 
solenoidal magnetic field. These devices detect charged 
particles and measure their momentum and energy loss. 
Photons and neutral hadrons are detected in a CsI(Tl) 
crystal electromagnetic calorimeter (EMC). The EMC 
detects photons with energies as low as 20 MeV and 
identifies electrons by their energy deposition. An inter- 



nally reflecting ring-imaging Cherenkov detector (DIRC) 
of quartz bars is dedicated to charged particle identifica- 
tion (PID). Penetrating muons and neutral hadrons are 
identified by the steel fiux return (IFR) , which is instru- 
mented with 18-19 layers of resistive plate chambers. 

The data correspond to a total integrated luminos- 
ity of 50.9 fb~^ taken on the T(4S') resonance and 6.3 
fb~^ taken off-resonance at an energy 0.04 GeV below 
the T{4:S) mass and below the threshold for BB produc- 
tion. In this sample, there are 55.5 ± 0.6 million BB 
events (Ngg). 

In this analysis, all charged track candidates are re- 
quired to have at least 12 DCH hits and transverse mo- 
mentum greater than 100 MeV/c. The track candi- 
dates not associated with a Kg decay must also orig- 
inate near the nominal beam spot. The muon, elec- 
tron, and kaon candidates must have a polar angle in 
radians of 0.3 < 9,, < 2.7, 0.410 < < 2.409, and 
0.45 < 6k < 2.50, respectively. In addition, all charged 
kaon candidates are required to have a laboratory mo- 
mentum greater than 250 MeV/c. These requirements 
ensure the selection of tracks in the regions where the 
acceptance is well understood by the PID systems. 

Photon candidates are identified from energy deposited 
in contiguous EMC crystals, summed together to form 
a cluster with total energy greater than 30 MeV and a 
shower shape consistent with that expected for electro- 
magnetic showers. 

Electron candidates are required to have a good match 
between the expected and measured energy loss (dE/dx) 
in the DCH, and between the expected and measured 
Cherenkov angle in the DIRC. The measurements of the 
ratio of EMC shower energy to DCH momentum, and 
the number of EMC crystals associated with the track 
candidate must be appropriate for an electron. 

Muons are selected based on the energy deposited in 
the EMC, the number and distribution of hits in the IFR, 
the match between the IFR hits and the extrapolation of 
the DCH track into the IFR, and the depth of penetration 
of the track into the IFR. 

Charged kaon and pion candidates are selected based 
on energy loss information from the SVT and DCH and 
the Cherenkov angle measured by the DIRC. 

The intermediate states in the indicated decay 
modes used in this analysis, J/ip {ee, fifi), 4> {K^ K^), 
77(77, Tr+TT^Tr*^), ry' (?7 (77) 7r"''7r~), 7r'^(77), and 
Kg(iT^Tr~), are selected with the mass intervals in 
Table U Since B° J/tpj] and B° -> J/^p7]' involve 
decays of a pseudoscalar meson into a vector and a 
pseudoscalar meson, the angular distribution is propor- 



ispr 

tional to sin^ 9i, where 9i is the helicity angle [ij of the 
lepton from the J/ip- Hence an additional requirement 
of |cos6'£| < 0.8 is applied to reject continuum and other 
backgrounds. The 77 candidates are rejected if either of 
the associated photons, in combination with any other 
photon in the event, forms a 77 mass within 20 MeV/c^ 
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of the 77° mass. For the mode J/tfjri{-f'y), the i] 

candidate is required to have |cos0^| < 0.8, where 0!^ 
is the photon hehcity angle in the rj rest frame. This 
rejects combinatoric background due to random pairs of 
photons that typically have a photon helicity angle that 
peaks at or 180 degrees. For the rj' ^ rj (77) tt+tt^ 
candidates, we use the same 77 selection criteria for the 
r] described above, including the 7r° veto. 

An additional requirement is applied to separate two- 
jet continuum events from the more spherical B meson 
decays. The angle Or between the thrust direction of 
the B meson candidate and the thrust direction of the 
remaining tracks in the event is calculated. We require 
\cos 9t\ < 0.8, since these thrust axes are uncorrelated 
and the distribution in cos^t is flat for BB events, while 
the distribution is peaked at cos0t = ±1 for continuum 
events. 

TABLE I: Mass regions for selection of intermediate particles. 

Mode Mass Range (GeV/c^) 

J/'tp e+e" 2.95 < M(e+e-) < 3.14 

J/ijj ^j.+ fi- 3.06 < M {fi+n') < 3.14 

<t)-^K+K- 1.004 < M{K+K-) < 1.034 

K%^Ti+Ti' 0.489 < M(7r+7r~) < 0.507 

77 ^ 77 0.529 < M (77) < 0.565 

77 ^ tt+tt-tt" 0.529 < M (71+71^71") < 0.565 

77' T^TT+TT" 0.938 < M (j77r+7r") < 0.978 

7r°^77 0.120 < M(77) < 0.150 



The intermediate candidates are combined to construct 
the B candidates for the six decay modes under study. 
The estimation of the signal and the background employs 
two kinematic variables: the energy difference Ai? be- 
tween the energy of the B candidate and the beam energy 
El in the T(4S') rest frame; and the energy-substituted 

mass mes = \J {El)^ — (Pb)'^, where Pg is the recon- 
structed momentum of the B candidate in the T(4S') 
frame. Typically these two weakly correlated variables 
form a two-dimensional Gaussian distribution for the B 
meson signal but not for background. The resolutions in 
AE and toes are decay mode dependent. A signal region 
for each mode is defined as a rectangular region in the 
AE versus TOes plane, hsted in Table HTl The TOes range 
is given in term of tties — tub-, where uib is the mass of 
B meson. The number of data events, no, observed in 
the signal region for each mode is listed in Table ^ 

The efficiencies for each mode are determined by Monte 
Carlo simulation. The simulations of J/ijj(j>K and J/ip(f> 
decays assumed three- and two-body phase space, respec- 
tively, with unpolarized J/ip and (j) decays. The J/ipV 
and J/tpT]' simulations used the angular correlations de- 
termined by the helicity amplitude. 

The backgrounds in the toes distribution are composed 
of two components: a combinatoric background, whose 



shape is described by an ARGUS function [7| , and a peak- 
ing background that peaks in the signal region and is 
described by a Gaussian function. The sources of com- 
binatoric background are the continuum events and two 
categories of BB events: decays with a leptonic J/4> de- 
cay, and those without. Monte Carlo simulation studies 
show that the source of the peaking background is BB 
events that contain a leptonic J/ip decay. 

The shape of the ARGUS function is determined mode 
by mode by fitting to the toes distribution of candidates 
in an enhanced fake J/ip sample, which is obtained by 
reversing the normal lepton identification requirements. 

The normalization of the combinatoric background for 
each mode is obtained from a fit to the toes distributions 
in the AE signal region of the on-peak data. The integral 
of the ARGUS function in the signal region is nc, the 
number of combinatoric background events. 

The peaking background is determined from a fit to 
the Toes distribution of Monte Carlo BB events with 
leptonic J/ ip decays using the sum of a Gaussian and an 
ARGUS function. The number of peaking background 
events np is the integral of the Gaussian function in the 
signal region. 

The total number of background events (n^) and the 
uncertainty on this number (ai,) are calculated from the 
fit value of riQ and np and their errors. The values of ni, 
and (Tfc are listed in Table^lfor all modes. The combina- 
toric background is by far the dominant background in 
all modes except the B'^ J/iI;i]{tt~^7t~tt'^) mode, where 
the peaking component is ~ 20% of the total background. 

In Table UTTl we list the contribution to the systematic 
error from the uncertainty on each of the following quan- 
tities: Ng-g; secondary branching fractions 0; Monte 
Carlo statistics; PID, tracking, and photon detection effi- 
ciencies, which are based on the study of control samples; 
and background parameterization, which is estimated us- 
ing AE sideband information. 

Additional systematic uncertainties due to the decay 
model dependence are estimated for the modes J/tp4': 
J/ip(f)K^ , and J/ip4>Kg. Monte Carlo simulations are 
used to determine how much the efficiency depends on 
assumptions about intermediate resonances and angu- 
lar distributions. Two samples are generated for each of 
the three modes with decay distributions determined by 
the assumed polarization of the vector daughter mesons, 
rather than by phase space. One sample is generated 
with 100% transversely polarized J/ip and mesons, and 
the other with 100% longitudinally polarized J/ip and 
mesons. The resulting relative change in efficiency is en- 
tered as a fractional systematic error in Table IIIII An 
additional check based on Monte Carlo samples with an 
intermediate state gives negligible effect. 

The total systematic error for each mode combines all 
these separate errors in quadrature and is listed (Total) 
in Table ITTIl 

There is evidence for signals in the J/^p(pK^ and 
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J /ip(j)K^ modes. The results are shown in Figs. |^ and 
13 The Poisson probability that the background nf, fluc- 
tuates up to the observed number of events, tlq, or higher 
is 7.7 X 10-6 for J/il)(t)K+ and 4.2 x IQ-^ for J/ipcjiK^g. 
The branching fraction for these modes is determined by 
a simple subtraction of events in the signal region that 
yields the number of signal events, Us = uq — ni,. The 
calculation of the branching fraction is based on the effi- 
and the secondary branching fractions 



ciency, Ug 



N 



for the J/i/j, 4>, and Kg from Ref. The results are 
summarized in Table Ull where the first error is the statis- 
tical error and the second is the systematic error, listed 
in Table HTH The derived result for B" J/iJjcIjK" is 
also shown in Table llll 
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FIG. 2: The AE and rriES distributions for B+ J/tp(l)K+. 
The AE vs. tties event distribution is shown in (a) with a 
small rectangle corresponding to the signal region selection 
defined in Table II. The AE projection with a mES signal 
region selection is shown in (b). The rriES projection with a 
AE signal region selection is shown in (c). The solid line in 
(c) is the fit described in the text. The Gaussian component 
includes both the signal and peaking background. 



For modes with no signal or limited statistical evi- 
dence (J/V'0, J/i^ri^ J/tpr]'), we determine both a cen- 
tral confidence interval and an upper limit interpretation 
for the branching fraction. The upper limit method uses 
no, Ub, and ab, in the signal region, and the total sys- 
tematic uncertainty ctt- Assuming the two uncertainties 
{ab, ot) are uncorrelated and Gaussian, the Bayesian up- 
per limit on the number of events (-/VgQ%) is obtained 
by folding the Poisson distribution with two normal dis- 
tributions for these two uncertainties and integrating it 
to the 90% confidence level (C.L.). In Table |ll| we list 
for each mode the efficiency, the number of observed 
events, the expected number of background events, the 
90% C.L. upper limit for observed events, the correspond- 
ing branching fraction limit and a central interval for 
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FIG. 3: The AE and mES distributions for ^ J/^p(j)K°s. 
The descriptions of Figs. 3(a), (b) and (c) follow those of 
Figs. 2(a), (b) and (c), respectively. 



the branching fraction. The upper limit obtained from 
the combination of the two B'^ — > J/^'V modes is shown 
in Table ^1 The mean value of the branching fraction 
is calculated for B° J H-'<P and J/''P'n'- We 

also combine the observed numbers of events for the two 
B^ — > J /ipr] modes to calculate a branching fraction of 
(1.6±0.6(stat.)±0.1(s2/si.)) x 10-^. The Poisson proba- 
bility that the background fluctuates up to the observed 
number of events or higher is 2.5 x 10^^ for the combined 
result. 

In summary, we determine the branching fraction of 
B J/ij(t)K in two modes, BiB+ J/ip(j)K+)= 
(4.4 ± 1.4 ±0.5) X lO"'"^ and B{B° J/V'0i^g)=(5.1 ± 
1.9 ±0.5) X 10"^ The branching fraction of B ^ J/^(j)K 
is consistent with a CLEO |lfl| resuh, (S.Slg;^ ± 1.3) x 
10^^. Upper limits have been determined for the modes 
B^ J/ip(l), J/iprj, and J/iprj'. The upper limit on 
B'^ —f J/ipr] is a significant improvement over the pre- 
vious best limit of < 1.2 x 10"^ at 90% C.L., from the L3 
Collaboration 0. The combined branching fraction for 
B^ J/4>r] is comparable to the B^ — > J/tpn^ branching 
fraction 1] . The search and resulting upper limits on the 
branching fractions for B^ J/4'v' and B^ — * J/4'4' are 
presented. 
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chine conditions provided by our PEP-11 colleagues, and 
for the substantial dedicated effort from the comput- 
ing organizations that support BABAR. The collabo- 
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TABLE 11: Branching fractions and 90% C.L. upper limits. 



Mode 


AE (MeV) 


Signal Region 

ImES-mg] (MeV/c2) 


Efficiency ng 


rib ± o"b 


-^^90% 


90% C.L. Upper Limit 
(10-5) 


Branching Fraction 
(10-5) 




57.0 


8.0 


10.6% 


23 


7.8 it 0.6 






4.4 ±1.4 ±0.5 




57.0 


8.0 


8.6% 


13 


3.3 ± 0.4 






5.1 ± 1.9 ± 0.5 


















10.2 ± 3.8 ± 1.0 




57.0 


8.0 


12.1% 


1 


0.3 ±0.2 


3.60 


< 0.9 


0.18 ±0.26 ±0.03 


J/H 


100.0 


10.0 


2.5% 





0.5 ± 0.3 


1.81 


< 6.3 


-1.7± 1.0±0.2 


J/ipr] (77) 


100.0 


10.0 


15.5% 


8 


1.7 ±0.4 


11.5 


< 2.9 




J/ipV {lT'^n~TT^) 


72.0 


10.0 


8.7% 


4 


1.5 ± 0.9 


6.76 


< 5.1 




J/ipV combined 














< 2.7 


1.6 ±0.6 ±0.1 



TABLE IIL Systematic error summary on the branching fractions. All are fractional uncertainties in percent. 



Mode 


Nbb 


Secondary 


Monte Carlo 


PID, Tracking, 


Background 


Model 


Total 




Branching Fractions 


Statistics 


Photon Detection 


Parameterization 






J/iP(l,K+ 


1.1 


2.2 


1.6 


8.2 


5.9 


0.4 


10.4 




1.1 


2.2 


2.1 


8.3 


1.9 


0.9 


9.3 


J/i"t> 


1.1 


2.2 


1.6 


6.7 


12.0 


1.0 


14.1 




1.1 


3.8 


4.6 


9.3 


7.1 




13.3 


J/ipv (77) 


1.1 


1.8 


1.6 


6.0 


6.9 




9.5 




1.1 


2.4 


2.2 


7.7 


8.0 




11.6 
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